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The r-function theory of Painleve systems is used to derive recurrences in the rank n of 
certain random matrix averages over U{n). These recurrences involve auxilary quantities 
CO ■ which satisfy discrete Painleve equations. The random matrix averages include cases which 

, can be interpreted as eigenvalue distributions at the hard edge and in the bulk of matrix 

^SJ \ ensembles with unitary symmetry. The recurrences are illustrated by computing the value 

' of a sequence of these distributions as n varies, and demonstrating convergence to the value 

'^^ appropriate limiting distribution. 

cn 



1 Introduction 



>■ ■ 1.1 Motivations and objectives 
O , 

C\J , In a recent series of papers ^S1~|IH1 we have shown how the Okamoto r-function theory of Painleve 

I systems can be applied to rederive known evaluations of certain random matrix averages in terms of 

' Painleve transcendents. Moreover it was shown how this theory could similarly be used to evaluate 

■ random matrix averages not known from previous studies, and to also yield recurrences of the discrete 

, Painleve type for the shift by unity of a parameter or parameters in the same random matrix averages. 

Q^l Subsequent to our works ^1 two different major theories — one on the discrete Riemann-Hilbert 

problem due to Borodin (51[7], and the other based on the integrable Toeplitz lattice due to Adler and van 
Moerbeke [5] — were applied in |SI El and |T] respectively to also provide recurrences for random matrix 
averages with respect to a shift by unity of a parameter. The averages considered were with respect to 
the unitary group U{n), and the shift performed in the rank n of the matrices. As with our own work, 
the average itself is related to an auxilary quantity or quantities, and it is the latter which satisfy the 
' primary coupled recurrences. 

' It is our objective in this work to further develop the Okamoto r-function theory as it relates to spec- 

ifying recurrences for random matrix averages. Whereas in our earlier works recurrences were obtained 
mostly with respect to an otherwise continuous parameter within the average, in the present work, as 
with the works by Borodin, and Adler and van Moerbeke, our attention will be focussed on obtaining 
recurrences with respect to the rank of the random matrix and thus the dimension of the average itself 
(the averages under consideration couple only to the eigenvalues of the matrix) . 

Typical of the results of this paper is the recurrence obtained in our work jl5| for the particular PIV 
r-function 



j = l l<j<k<n 
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/ GUE 
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Here GUE„ refers to the probability density function 

n 

cH^"'^' n (1-2) 

j=l l<j<k<n 

with C denoting the normalization, realized by the eigenvalues of Hermitian matrices with certain 
complex Gaussian entries (see e.g. ^T]) and Xj"* = 1 for Xj G J, x'P = otherwise. From ^1 
eqs. (2. 86), (2. 4), (2. 75)] we have that for an appropriate C„ independent of t (according to H2.11|l be- 
low Cn = 2n), and with T'^[n] T^^[n]{t: //), 

^" , Jr n ^ = 2n+ (2i - /o[n] - Mn])f2[n], (1.3) 

where /o [n] , /2 [n] satisfy the coupled recurrences 

2n 

fo[n] + Mn-l] = 2t~f2[n] + —^, n = l,2,... (1.4) 

/2[n+l] + /2[n] = 2t-/o[n]+ ^^" + ^ + ^\ n = 0,l,... (1.5) 

These coupled recurrences were shown to be equivalent to a single second order difference equation known 
in the literature as the discrete Painleve I equation. Specification of /o[0], /2[0], t^^[0], t^^[1] (sec (|2.15|) 
below) uniquely determines {/o['^]}n=i,2,..., {/2M}n=i,2,... and {T^^[n]}„=2,3,...- As noted in Section 2 
below, the more general PIV r-function 

n 

r^^N(t; = ( 11(1 - ex|lL))(i - ^^r)^^^^^ (1-6) 
j=i 

also satisfies the system H1.3|l ~ ()1.5|l . 

Recurrences with respect to the dimension of the random matrix will also be given for three averages 
over the unitary group U{N), known from our earlier work to be r- functions for certain Painleve systems. 
With zi := e**' these are 

N 

/uiN) 



mt;f,) := (n<^^^^"^^' %(^) (1-7) 

N 

r^[7V](t;M,i') := + ^^'Ti^ + ^ ^^iTe'^X.^. (1-8) 

1=1 ' ' 



1=1 

where U{N) refers to the probability density function 



n \^k~zA', (-vr<0, <7r, j-l,...,iV). (1.10) 

In the case of (|1.7f) we only take the U{N) average as the definition for /i G Z; for general ^ the r-function 
r^ii [N]{t; fi) is to be defined as the Toeplitz determinant given in 1)3.21) below. Also, as written H1.8|l is 
only well defined for /i, G Z>o. However with z = e'^, use of the identity 

(1 -I- Z)^(l -I- 1/z)'' = zt^-'^)/^!^ _|_ (1 11) 
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gives 

n 

T''[N]{t;f,,,.) -.^ {l[zl^"'^/'\l + zir+'^e'^') (1.12) 



UiN) 



which is well defined for Re(/i + i/) > —1. 
We also indicate how the PV r-function 



f^[n](i;/i,a;0 '-^ {Ri^ ' ^x\Z)){x, ~ tr) ^^^^ (1.13) 
j=i 



and the PVl r-function [HI 

n 

f^^N(t;^,a,&;0:=(n(l-^xSl\))(i-x,r)j^^^ (1-14) 
can be characterized by recurrences. Here LUE„ refers to the probability density function 



1 " 

.7 = 1 



j = l l<j<k<7i 

while JUE„ refers to the probability density function 



n 

j^nxgi)-."^!--^-)'' n (-^ (1-16) 



i=l l<j<k<n 

The normalizations in p.l5|l and are 

N N-1 



/•OG /'OG l y —± 

lN{a):= dxi--- dxNY[x1e-^' [] {xk - x^f ^ H T{2 + j)T{a + 1 + j) (1.17) 

■^0 -'0 J — 1 n 



l<j<fe<7V j=0 



and 

JN{a,b) := 



/ dxiXi{l- Xi)'' ■ ■ ■ dxN X%{1 - Xn)'' Yl i^k - Xj) 



l<j<k<N 

^^ r(a + i + j)r(b + i + j)r(2 + j) 

- n T{a + b+l + N + j) • ^^-'^^ 

We remark that both p.l5|l and (|1.16|) can be realized as the eigenvalue probability density function for 
certain ensembles of random matrices (see e.g. jll)'). We have not been able to derive recurrences for 
and (|1.14|) in n only; rather the recurrences to be indicated also act on the parameter a. 

1.2 Strategy 

The Okamoto theory is based on a Hamiltonian formulation of the Painleve equations, which in turn 
can be traced back to Malmquist [22]. Corresponding to each of the Painleve equations PII-PVI is a 
Hamiltonian if, which is itself a function of the conjugate variables p and q, the independent variable t, 
and a number of parameters. The conjugate variables p and q are also dependent on the independent 
variable t and the parameters. By eliminating p in the Hamilton equations 

, dH , dH ,^ 
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the Painleve equation in q results, although we have no explicit use for this defining feature of H below. 
A particular parameter n is distinguished and w6 write H — , p — p^i , q — q^i . Our primary concern 
is in so called Schlesinger transformations, which relate the Hamiltonian system with parameter n + 1 to 
the Hamiltonian system with parameter n. 

One introduces a r-function t„, a function of the independent variable t and the parameters, by the 
requirement that 

ff„ = ^logr„. (1.20) 

^From the Okamoto theory it is known that 

^-^^f^=fiPn,qn,t) (1.21) 

for some explicit polynomial function /, typically related to the time derivative of Furthermore, the 
Schlesinger transformation theory gives that {p„, qn} satisfy coupled first order recurrences 

Pn+l = 9l{Pn,qn), g„+l = ^2 (Pn , ?«) (1-22) 

for some explicit rational functions 51, (?2- Thus once po,qo have been specified {pm Qn}n=i,2,... can be 
generated from (|1.22|) . With this information, and knowledge of to,ti, H1.21|I can be iterated to specify 

{t2,T3,...}. 



1.3 Paper outline 

We will devote separate sections to each of the r-functions (|1.6|I - H1.9|I . with (|1.13|) and H1.14|l considered 
during the discussion of (|1.8|l and H1.9|l respectively. In the cases of (|1.6f) . (|1.7|) and (|1.8|) the Schlesinger 
transformations which increment the dimension of respective random matrix averages are known from 
our earlier works |15l I16|. The formulation of the recurrences is then a straightforward application of 
the strategy outlined above. However in the case of H1.9(l there is some complication as one must first 
change variables to obtain a random matrix average for which the standard Schlesinger transformation 
increments the dimension of the random matrix average. In the final section some uses of our recurrences 
for the computation of the random matrix averages as they occur in applied problems will be discussed. 



2 The r-function sequence {r^^[n] /i, ^)}„^o,i,... 

The Hamiltonian for the PIV system is given by [221 

H'^ = {2p ~ q ~ 2t)pq - 2aip ~ a2q. (2.1) 

Let 

{ai,a2)^{~H,-n) (2.2) 

and write H^^ = H^^^ thus distinguishing the parameter a2 = —n. It was shown in ^1 Prop. 22] that cor- 
responding to the sequence of Hamiltonians {-ff^^}n=o,i,... is the sequence of r-functions {r^^[n](i; n)}n=o,i, 
as specified by (|1.1|) . Moreover, combining the result of ^1 Prop. 6] with the workings leading to ^1 
Prop. 7 and Prop, 22] it follows that more generally T^^[n]{t; iJ,,£,) is a r-function for iJ^^. 
The significance of this latter fact is that with 

fo[n] 2t + g„ - 2p„, /2M 2p„ (2.3) 
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we know from 15, eq. (2.75)] that the recurrences (|1.3|l - H1.5ll hold, and these recurrences fully determine 
{T'^[n]{t;fi,^)}n=2,3,... once we specify /o[0],/i[0] in and C„, r^^[0], t^^[1] in (O. To 

determine C„ we require the fact ^| eqs. (2.41), (2.42)] that with 

C„ = ^"+^^"-\ jr.e'W[n]^a'^H r^^[0] = a^^[0] = 1 (2.4) 
the function has the explicit double Wronskian form 

a^^[n]=det|^a^^[l]l , . (2.5) 



jM=0,...,n-l 



Now it follows from jl5l eq. (2.41), Prop. 6] that up to a proportionality constant, which we are free 
to choose to be unity, 

. noo POO 

a^^[l] = e*'(/ ~U )it-x)^e-''' dx. (2.6) 



Noting that can be written 

/•OO POO 

a'^[\] =(^J -ij y~xre-^"-^'^dx 
the differentiation required by (|2.5() becomes simple to perform and we obtain 



dt'+i 



OO pOO 



cr^^[l] = 2'+Je*'( / yt - x)''+'+^ e-""^ dx. (2.7) 



Substituting this in H2.5|l and recaUing the workings of ^1 proof of Prop. 21] we see that 

poo POO POO POO ^ 

W )^^i-"(y -^J^ )dx,,l[e-^Ht-x,Y 



nn(n — l) 

ni 



OO 

X n {xk-xj)^ 

l<j<k<n 



It is well known (see e.g. ^J) that the normalization C in the definition p.2|l of the GUE„ probability 
density function has the explicit form 



C = n!2-("-i)"/V"/2 J] U 

1=0 

so (I2.8|l can be written 

n — 1 n 

;=o 1=1 

n-l 

= 2"("-i)/2^"/2 J]/!e*'V^[n](t;^;0- (2.9) 

1=0 

Recalling H2.4|) we thus have 

n-l 

7„ = 2"("-i)/27r"/2 Jl (2.10) 
/=o 

and this in turn implies 

Cn = 2n. (2.11) 
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Regarding the initial conditions for ()1.4|l and (|1.5(l . we require the facts |[15l proof of Prop. 6] that 

Po = 0, go = ^logT^^[l]. 
Thus recalhng \2.'S\\ and H1.6|l we have 

/o[0]=2i+-log(( / / ){t-xYe-^' dx), /2[0] = 0. (2.12) 

The initial conditions for H1.3|) are by definition 

t'^[0] = 1, r^^[l] = -^(y -ey^ )(t-xre--'dx. (2.13) 

In summary, we thus have that the foUowing resuh holds. 

Proposition 1. Let r^^[n] = t^^ [n](f; /i; ^) as specified by 1^1. b]) . Let p^, (?„ denote the conjugate 
variables in the Hamiltonian \2.1]) with parameters given by ^2.1^} . and define fo[n\ and in terms of 
these variables by 1^2. We have that {/oM}n=i,2,..., {/2M}n=i,2,... md {T^^[n]}„=2,3,... are determined 
by the recurrences J^1.3\) - J^l.^) subject to the initial conditions ^2.12^) . Y2.1S\) . 

We remark that in the special case ^ = 0, /i G Z>o, Hl.till is a polynomial in t, which in view of (|2.4|l . 
1)2. 5|l . (|2.7() and the integral representation 

{t - ix)Pe-'=' dx = V^2-PHp{t) 



(20~^r(^+'-)i/^+j+fc(it) . (2.14) 

-I j,k—0,...,n—l 



has the explicit form 

T^^[n] = — det 

In 



In this case ()2.12f) and H2.13|l can be written 



r^^[0] = l, r'V[l]^^^H,i^t). (2.15) 

71 

It is of interest to recall the duality formula ^1 eq. (4.37)] 

r^^[n](t;A*,0) = z-"^T^^M(zt;n,0) = — det[2-"z-(^+'=)H„+,+fc(i)],- fe=o,...,M~i (2-16) 

thus giving (|2.14|) for n = 0, 1, . . . as a sequence of /i x /i determinants. 

Another point of interest is that with the initial conditions (12.12(1 a closed form solution of the coupled 
recurrences (|1.4f) and 11.5|l can be given. Thus it follows from 24, eq. (4.8)] that 

f r . y[n]it:ti,OT'''[n + l]it;i,+ l,0 r^^[n + l](t; /i, Qr^^jn - l](t; ^ + 1, Q 

TIy[n]it■,^, + l,Or'''[n+m■,^i,0' ^'^""^ r'y[n]it; Or'^[n]{t; + 1,0 

(2.17) 

(the proportionality constants cannot be read off from |24| : these are determined by considering the 
t 00 behaviour of (|1.3|l - p.5|l '). 
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3 The r-function sequence {r^-^'^'[A^](t; /i)}Ar^o,i,... 

Although (|1.7|l is well defined for all complex ^, only for /i e Z will we take the U{N) average as the 
definition. For general fi we will make use of a Toeplitz determinant form, obtained by applying the well 
known identity 



N 



1=1 



This gives 



N 

n 

1=1 



U{N} 



U{N) 



det 



2tt 



j,k=l,...,N 



(3.1) 



det 



r 1 

2^ 



g«(M+j-fc)9gVtcose ^0 



j,k=l,...,N 



,N- 



(3.2) 



The second equality of H3.2|l follows from an integral formula for I,^{z), valid for i/ e Z. We take this 
latter determinant as the meaning of r^^^ [N] for general /j,. 
Now the Hamiltonian for the PHI' system is given by |2Hj 



tH 



III' 



q p ~ {q + viq - t)p + -{vi + V2)q. 



We showed in 16, Section 4.3] that with 

{vi,v2) = i^I + N,-^l + N) 

the quantity 



r''^/^ det[i^+,^k{Vt)U=i., 



(3.3) 



(3.4) 



(3.5) 



is a r-function for the corresponding sequence of Hamiltonians {tHjJ^ }7v=o,i.... (it is still valid to call 
(|1.7|l a r-function for a PHI' system as the extra factor is equivalent to the addition of a constant 

to the Hamiltonian H3.3() . which of course does not alter the Hamilton equations). From the working in 
[T7] summarizing the Okamoto theory of PHI', we can deduce the following recurrences for (|1.7|l . 

Proposition 2. LetT'"'[N] = t'"' [N]{t; ^i) as given by jl.T}) , and let pN, qjv denote the conjugate 
variables in the Hamiltonian f!^.yp with parameters j^^p. The sequences {t^^^ [-^]}a/'=o,i,...; {Pn}n=o,i,...! 
{9^1^=0, 1,... satisfy the coupled recurrences 

t"^' [N + iy^^' [N-l] 



PN+l 



iJ-qN , 



PN 



1 



t 

qN+1 — 



{1 + N)t 



qN qNiqNipN - I) - tJ-) + t 



(7V= 1,2,...) 

(iV = 0,l,...) 

(TV = 0,1, 



subject to the initial conditions 



r^^O] 



qo^t-\ogt-^^H^{2Vt) 
'[l]=/^(2Vt). 



(3.6) 

(3.7) 
(3.8) 

(3.9) 
(3.10) 



Proof. The working of [161 proof of Prop. 4.2], which in turn is based on ^311, together with ^1 
eqs. (4.13), (4.20), (4.21)] tell us that 



^iii' 



[N + 1} 



-III' 



N - 



d_ 

di 
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so (I3.6|l now follows from H3.3|l . Furthermore, it was shown in |16[ eqs. (4. 40), (4. 41)] that for the Hamil- 
tonian 1)3. 3|l with parameters (vi,V2) — + n,v^2''' + 

Pn+l = f{Pn-l)~'fy^^-vi'^) + l (3.11) 

Qn+i + — A . 3.12 

Setting vf'^ = —v^2^'' ~ fi gives (13 .71) and (|3.8I) . The initial conditions (|3.9() follow from ^1 proof of 
Prop. 4.3], while the initial conditions H3.10|) are immediate from H3.2|) . □ 

It is known (see ^1 Prop. 4.6]) that the sequence {qN}N=o,i.... satisfies the particular (alternate) 
discrete Painleve II equation 

^ + ^ + ^ ^±_q.^N-,^ ^^0,1,... (3.13) 

gNQN+l+t qN-lQN + t Qn t t 

In the special case /i = the sequence {pn}n=a.i,... itself is also determined by a particular discrete 
Painleve II equation. To see this, note that 1)3. 7|l with /i = gives 

2 , 1 ^ PN+i , .^ 

qpf^t- . (3.14) 

1-pN 



Setting 

n\T = \/f.- 

Pn 



IN = VI^, Po = 1 (3.15) 



we see that (|3.14ll has the unique solution 



PN^l- Pn- (3.16) 



Making use of this in Proposition [3 we obtain the following recurrence scheme for {t'" [N]{t,0)}, first 
derived by Borodin 5 from a discrete Riemann-Hilbert approach, and subsequently obtained by Adler 
and van Moerbeke [J from their theory of the Toeplitz lattice, and by Baik P who used a Riemann- 
Hilbert approach distinct from that of Borodin. 



Proposition 3. We have 



1 - = 



T 



^^^'m + i]t'"'[n -I] 



where {Pn}n=i,2,... satisfies the particular discrete Painleve II equation 

Pn 4- 1 4^ Prt — 1 — — — , 71 = 1 , 2 , . . . 

^ ^/^(l-P2)' ' ' 

subject to the initial conditions 

Po^l, Pr = ^^. 

/o(2^/^) 

Let us show how g^r, like pn, can be written in terms of the r- functions. Put 

(t-^'^/V"^'[7V](i;/i))|t^4t =:r^ 

and denote the corresponding Hamitonian by H'^. Denote by Ti (T2) the Schlesinger operators with the 
action on the parameters (wi, W2) 1^ (wi + 1, t'2 + 1) ((wi , W2) (wi + 1, i'2 — 1)). It is known [2H| that 

T,tH^ = ti?^^, = tH^ + q^il - pn), T^tH^ = tH^+' = tH% - q^Pn (3.17) 

where q^ q^, PN ■— P^- Using 1)3.17(1 together with (|3.11(l . H3.12|l the sought formula can be deduced. 
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Proposition 4. We have 



log (%?^) = -tiH^nXl Ky + m = 



Proof. We have 

" - / ' n+l 'n 
'n+l ^" 

According to l|TTT|) . (jXHIl 



Thus 



qN+iPN+i - qNPN- 



qN+iPN+i = -qNiPN -l) + vi + l . 

qN 



= -^(2q%PN-q%-{vi + l)qN+t) 

™ ^ Tit ^ qN ^ ' 

1 fdtH'"' \ d , . , ^ 

= qN] = -t—\og{qN/t) 

qN ^ opN ' at 

where to obtain the final equahty, use has been made of the first of the Hamihon equations. This imphes 
(|3.18|l up to a proportionaHty constant. To determine the proportionahty, say, we use the asymptotic 
formula ^1 proof of Cor. 4.5] 

det[/,_fe+4^/^)],,,.0,...,„-l - ^nVi-inVmost+oii) 

t — ^oo 

which in light of H3.2|l and H3.6|l implies 

PN ^ (4i)-i/2 {N^Q) (3.19) 

t — ^oo 

while H3.18|l (with the proportionality still unknown) implies 

qN ^ cnVI. (3.20) 

t — *oo 

Substituting in H3.8|l and taking into consideration l|3.9(l implies cat = (—1)^. □ 



4 The r-function sequence {t [N]{t; fi,iy)} 



The definition H1.12|l of T^[N]{t; fi, i/) is well defined for Re(/i + ly) > —1. This domain can be extended 
by using H3.1|l to rewrite p.l2|l as a Toeplitz determinant and evaluating the integral, 



r^[7V](t;/i,i.) 



det 
det 



In 



r(^ + i/ + i) 



r(^ + j-fc + i)r(j/-j + fc + i) 



j,fc=l,...,7V 

iFi (-1/ + j - A:, ^ + 1 + j - k; -i) 



j.k=\.,...,N 

(4.1) 



Here the integral evaluation, which is well defined for general complex /u, follows by expanding the 
exponential in the first Toeplitz determinant and evaluating the resulting integrals using the formula 



1 

2^ 



r(a + b+ 1) 

r(a + i)r(6 + i)' 



(4.2) 
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In ^1 display eq. below proof of Prop. 3.6], {t^ [N]{t; fi, v)} n=o,i,... has been identified as the r- 
function sequence corresponding to a particular Schlesinger operator for the PV system. However some 
technical details of the derivation given there leads to complication for the present purposes, which in 
fact can be avoided by revising some of the workings in Let us then undertake such a program. 

We will construct a r-function sequence relating to the Hamiltonian |32[ 1^ 

iiJ^ ^ q{q - l)p{p + t) - {v2 - vi + V3 - V4)qp + {v2 - vi)p + {vi ~ V3)tq (4.3) 

for which eliminating p in the Hamilton equations shows 1 + l/{q — 1) satisfies the PV equation. Since 
eliminating p does not give the Painleve equation in q itself, we refer to this as the PV* system. (In 
we made use of the mapping between the PV and PV* systems |32) . which is in fact unnecessary and 
is what leads to complications for the present purposes.) Our interest is in the particular Schlesinger 
transformation with action on the parameters 

1113 

TQ^ivi,V2,V3,V4^) = {Vi - -,W2 - -,V3 - -,^4 + -). (4.4) 



We know from ^3 eq. (2.16)] that 



(4.5) 



This motivates introducing the sequence of Hamiltonians 



ttV ttV 



v^{v{°> -n/4,t)f> -n/4,D<°' -n/4,i;f' +3n/4) 



and the corresponding sequence of r-functions r„ specified so that 

I log (4.6) 

Following '77', the seed initializing the sequence of r-functions is a classical solution to the PV* 
system. 

Proposition 5. Let v^"^ — v^f"^ = 0. Then the PV* system admits the solution 

= (.(")- 4°) )t, ro = e(^i"'-^"')*, po.i|logrr + (4°^-4°))i (4.7) 

where e^^^i '^^3 '"'^tY satisfies the confluent hypergeometric differential equation 

ty" + iv<f> - + 1 + t)y' + (vf^ - v^^^)y = 0. (4.8) 

Proof. Direct substitution oi qo = I, v^'' — v^^ = into (|4.3I) gives the stated value of tH^ . The final 
equation in (|4.7|l follows from (|4.5|) . H4.3|l and H4.6|l which together give 

T,-hHr tHf ^ i| logrr* = tHf + Po ^ (4°^ - «f )i +Po. 

Now that the final equation in (|4.7|l is established, H4.8|l can be derived from the second of the Hamilton 
equations (|1.19|) 

dtH^ 



tp'o = = -(Mpo + t)- (.f - + (^"^ - (4.9) 

"3 "4 

by substituting the former equation for po throughout. □ 



10 



According to proof of Prop. 2.2], with 



(4.10) 



and Tq as in 1)4. 7|l so that tq = 1 (recaU that in H4.7|l we require — U4 = 0), the sequence {t„}„=2,3, 
is specified by the determinant formula 



T„ = dct[S' + ''Ti]j^k=0,...,n-l- 



(4.11) 



For an appropriate choice of the soKition of (|4.8|l and thus of ti , (|4.11() can be related to the r-function 
sequence (|1.8ll . 

Proposition 6. Of the two linearly independent solutions to i4.f^ , choose the solution analytic at the 
origin, 

(4.12) 
(4.13) 

(4.14) 



with 



JO) _ JO) 



(0) (0) 



Then i4.1U\ ) and give 



^'^tv'_( r(/i + i) v 



r(/i + i/ + i) 



Proof. Choosing tY as in (|4.12|) and the parameters as in 14.1311 gives, upon comparing with H4.1|l 



, y, r(/i + i)r(z. + i) _^ 



r(^ + i/ + i) 



Furthermore, use of H4.10|l and (|4.11|) shows 



'r(/i + i)r(j^ + i)\» 



V 27rr(^ + 1) 
We know jlHL proof of Prop. 3.1] that H4.11|l is equivalent to 

i"^f„ = det[,5J'+'^-i^fi],- fc^o....,„-i 



0,...,n— 1 



(4.15) 



(4.16) 



for any c, and thus choosing c = —1/2 the factor of t^/'^ in the determinant is cancelled, while the left 
hand side of H4.15)l is multiplied by t^"/^. Now simple manipulation involving integration by parts shows 



(5(e*y (l + z)^(l + l/z)'^e*^d6') 

= e'({v + l) [ {l + zY{l + l/zY+^e''de-{ii + y+l) f (1 + z)^(l + 1/z) V" d^) , (4.17) 



so making use of 1)1.11(1 and applying H4.17|l to column k and subtracting (/i + + 1) times column fc — 1 
for fc = ri — 1, 71 — 2, . . . , 1 in order shows H4.15|l can be reduced to 



X det 



6^e' 



{1 + zY{l + l/zye'' d9 (5^+'^-ie* / {1 + z)^'{l + l/zY+'e'' dd 



j=0,..-,Tl-l 

k = l,...,n-l 



11 



Repeating this procedure for columns fc = n— l,n— 2,...,2 and so on shows 



n(n-l)/2 (v+n)t V* 



r(A^ + i) 



n-1 



27rr(^ + i/ + l) 



-j \{V{v + l + l)dei\5'e' [1 + zY [1 + I / zY+*' e'^ dO 



1=0 



j,k=0,...,n-l 



and apphcation of the general identities [TH] 



det 



det 



j,fc=0,...,n-l 



j,fc=0,...,n-l 



= (u(i))"det[<5Vfe(t)],,fc=o.....n- 



i,A:=0,...,n-l 



then gives 



0,...,n — 1 



(4.18) 



The stated result now follows after noting that the factor (1 + l/z)*^ in the integral can be replaced by 
(l/zY without changing the value of the determinant. □ 

Knowledge of (|4.14|) and key recurrences from the Okamoto theory of PV as detailed in TB' allows 
the following recurrence for T^[N]{t; ^, v) to be deduced. 

Proposition 7. Let t^[7V] = T^[7V]{i; n, v) as given by il.^) or more generally j^. Letpjs[,q]sr denote 
the conjugate variables in the Hamiltonian i4.S{ ) with parameters 



vi -vs = -ly, vi - V2 = H; V4 -V3 = N 



and define 

The sequences {T^[iV]}Ar=o,i,...; {xn}n=o,i,..., {yN}N=o,i,... satisfy the coupled recurrences 

.T^[N +1]t^[N -1] ( t M\ / 1 

{T^[N\y \ yN 2/\yN 



1. X 1 

XN (Pn + t)qN ~ -;;(V2- vi), yN ^ — ■ 

2 qN 



t 



xn + xn-i = 

I 

yNyN+i = t 



VN 

N 



N 



yN 1 - yN 
XN + iy + ^i/2 + N + l 

x% - (m/2)^ 



subject to the initial conditions 



xo + fi/2 + t—\ogiFi{-iy,n+l:-t), yo = 1, 
at 



[0] - 1, 



r(/i + i)r(i/ + i)- 



Proof. According to ^1 proof of Prop. 2.2] 



V' v 
'n+l'^'n-l 



irX'r dt 



(4.19) 
(4.20) 

(4.21) 
(4.22) 
(4.23) 

(4.24) 
(4.25) 

(4.26) 
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so taking into consideration (|4.14() . (|4.3() and (|4.19|) we arrive at H4.22|l . H4.23|l . Furthermore in ^ 
Prop. 2.4] it was shown that with xn^Vn specified in terms of pn, Qn by (|4.2U|) . {xnjUn} satisfy the 
discrete Painleve IV recurrences 



VnVn+i = t 



t w 



(0) „(0) 



(4.27) 



Making use of (|4.19(l then gives 14.22|l and 14.23(1 . The initial conditions follow from H4.2()|l . 14.7|l and 

KWi . a 

For n e Z>o, the confluent hypergeometric fimction iFi(—n;c;~t) is proportional to a Laguerre 
polynomial, 

, r(a + l)r(n + 1) , 
T(n + a + 1) 

Thus for V e Z>o it follows from (lO) that 

r^[iV](t;/i,z.) = detfcr-M)l , „ 
(note that L^{-t) for n < 0). According to (tO^ we also have 



(4.28) 



a;o = t + + /i/2 - (i^ + /u)- 



(4.29) 



As in the PIV theory, in the case ly G Z>o the N x N determinant for T^[iV] can also be expressed 
as a v X v determinant. Thus we know from 1 16!, Props. 3.6,3.7] that for G Z>o 



r 

r^[7V](i;A*,i^)cxdet ^_L^^,(-t) 



Using the Laguerre polynomial identities 

Lri(x)=L^(x)-L^i(x), 

this is equivalent to 



j,k=0,...,u-l 



(4.30) 



dx 



L^ix) = -L;^t{x) 



'[N]{t;fj,,iy) oc det 



L 



fj.+j-k 
N+k-j 



i-t) 



j,k=l,...,ii 



and thus 



'[N]{t■^l,v)^r^[v]{t■^l,N). 



To determine the proportionality, we use the fact, following from p.8|l and that 



where 

Mjv(a,6) 



"'-1/2 J -111 t~; i^.tt^« 



'-1/ 

n 



1/2 

r(a + b + i + j)r(2 + j) 
r(a + i + j)r(5 + i + jr 



l<j<k<N 



(4.31) 
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Hence for v e Z>o 

''l"l('--'^M^''H«:.,A'). (4.32) 

Another point of interest is the exphcit r-function form of the sequences {a;Ar}Af=o,i....i {j/Af}Af=o,i.... 
generated by the discrete Painleve IV coupled recurrences (|4.22l) . (|4.23|) with initial conditions 14.2411 . 
We find 

/o + M + 1' ^)^[^](^; /i - 1, + 1) 



T 



[Ar + l](t;M, r.)T[^](t;M,J^+l) 



r[N]{t■,^,~l,:.+ l)T[N]{t■,^i + ^ 

These formulas can be established by proceeding in an analogous fashion to the proof of Proposition 01 
Let us now turn our attention to the PV r-function as specified by In some special cases 

this is intimately related to as specified by H1.12(l . Thus with In (a) specified by (|1.17(l we have ^1 
Prop. 3.7] 

-T [N]{-t;n,a;0) = —— — — r [^](t;a,7V), E Z>q (4.34) 



iNia) ^ ^^f'^}y [a]{t■,^,,N), aeZ>o (4.35) 



The identities H4.34|l and H4.35|l allow those special cases of f ^ [n] to be computed by the recurrences 
of Proposition 13 however the recurrences will no longer be with respect to the dimension of the average 
n but rather with respect to one of the parameters. For general parameters a system of recurrences for 
'r^[n] can also be given, but these recurrences alter both n and the parameter a. This comes about as a 
consequence of the following analogue of Proposition El 

Proposition 8. Write the general solution of the confluent hyper geometric equation \4-<^ in the integral 
form 



and set 



We have 



(0) (0) , 1 (0) (0) 



n-l 

^ „(7-l)t+-»(7 



i)(r(7~a))"-i nr(^ + i)r^N(i;«-i,7-a-^;C)- (4.36) 



i=0 

Proof. With 



F{a,'y;t) -.^e'i^J J ) e-*"(u - l)"-iu'^~"-i 
we see from H4.1()|l . ()4.11(l and 14.16|l that 

Analogous to (|4.17|) we can show that 

SF{a, 7; t) = -aF{a, 7; i) - (7 - a - l)F{a + 1, 7; 0- 
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Then proceeding as in the derivation of H4.18|l we deduce 

e-^-^-DVr ^ ~ det[e-*f (« + fc, 7 + j; t)kfe^o,...,n-i. 

n;=i r(7-a-0 

The method of the proof of ^3 proof of Prop. 3.1] allows this determinant to be written as a multiple 
integral. Making use too of (|1.17|l then gives 14.36() . □ 

Because (|4.26f) and (|4.27|l hold for any r- function sequence with the property 14.4|l . we see that (|4.36|) 
allows us to specify the analogue of Proposition [3 for {f^[n]{t; a — 1, 7 — a — n; ^)}n=o,i,..., although we 
stop short of writing it down. 

5 The r-function sequence {T^^[n\{t; ii,wi,W2; ^)}n^o,i,... 

As written H1.9|l requires — tt < ^ < tt to make sense, however we can readily extend this definition to 
general complex 

t := e*^ (5.1) 

First we make use of (|3.1|) to obtain 

T^^[iV]=det[i- /^l-ex(.-0..))e™='|l + ^P"He^^""^^-^P^^^-'rf^l ., , (5.2) 



j,k=l,...,N 



The integral in H5.2|l naturally breaks into two. Introducing t according to ()5.1|l and setting n £ Z>o the 
first portion reads 

where, with w := wi + iw2-, 

a = iv~^, + j — k, b = w + ^ — j + k 

(the integral evaluation follows upon making use of (|4.2() '). For the second portion, taking both wi,fj, € 
Z>o, and writing in terms of dz instead of d6 we have 



27ri 



•'C(_i/t__i) ^ 



where C(_i/f._i) is a simple closed contour starting at z = —1/t and finishing at z = —1. Making the 
successive transformations z s- —z/t, z h- s- — z +1, z 1-^ (1 — t)z, then making use of the integral formula 

, ' x^^ (1 - x)^-^ (1 - txy dx = ^^t]^^}^}^Yt,^^ 2F1 (r, Ai + 1, Ai + A2 + 2; t) 
shows (|5.4ll is equal to 

C^~^^ ±7ri(fc-7+Ai-tB) ^(^A^ + l)r(2wi + l) _^fc--,-|-^— ffi^ _ ^\2p+2!iii+l 

27ri r(2^ + 2wi + 2) ^ 

X2F1 (2^ + 1, 1 + fc - j + ^ + ly; 2^ + 2wi + 2; 1 - i) (5.5) 

where the ± sign is taken accordingly as Im(t) 0. Substituting for the integral in H5.2|l the sum of the 
hypergeometric functions (|5.3|l and (|5.5|l gives meaning to T^^[Af] for general complex values of H5.1|l . 
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We know from that the CUE at average 1)1. 9|l can be written as an average over the generaUzed 
Cauchy unitary ensemble [331 E| specified by the p.d.f. 

1 ^ 1 

where Afjv is given by H4.3f |l . Thus making the change of variables 

1 + ixi 

zi — 

1 - ixi 

in shows [TTI eq. (1.19) with ^ i-^ 2/i] 

(5.7) 

where 



r 1 - (1 - C)e~"'^ (5.8) 

Furthermore, we know from QX^ that the CylJE average in H5.7|l as a function of s is also the r-function 
for a particular PVI system. Moreover, unlike the situation with (|1.9|l . the Schlesinger transformation 
studied in 1 7 increments N in this average and leaves the other parameters unchanged (in (|1.9|) this 
same Schlesinger transformation increments N but also decrements fi). 

To make these statements more explicit, we recall |2E] that the Hamiltonian for the PVI system is 
given by 

t{t - l)i7^^ = q{q - l){q - t)p^ - ({vi + V2){q ~ l){q - t) 

+ {vi - V2)q{q -t) + (W3 + Vi)q{q - 1))]3 + (wi + Vi){vi + Vi){q - t). (5.9) 
Introduce the Schlesinger operator T3 with action on the parameters 

J3V = {vi,V2,Vz + 1,^4), 

and with appropriate actions on the conjugate variables p, q. We know [171 eq. (2.27)] that 



rpriTTVI . TjVI TjVI 



/ (0) (0) (0) , (0), ■ 



^From this we introduce a sequence of r- functions t„ specified so that 
We know from 17", Prop. 15] that 

N 

CyUE„ 



1=1 

where with jj := rji + ir]2 



n(i-«*!lU)(--)»>,,„,J„_„,„-S'(T;v) ,5.10) 



V = (-771, i772, 771 + iV, + ryi). (5.11) 

For our present purposes the proportionality constant in (|5.10() . not calculated in |17|. is of importance. 
To obtain its value we must recall some of the results from jl7| . 
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We know from (ni eqs. (2.30), (2.37)] that with 

fy := {t{t - i))(«+-r +.^'")(n+4"'+-r)/2^v/ (5.12) 

we have 

f„^^ = det[5^'+'=frkfc=o.....«-i, 5:=t{t-l)^. (5.13) 



We also know that t^^ satisfies the Gauss hypergeometric differential equation 

t{\ - t)y" + (^c - (a + 5 + l)t^y' ~ aby ^ (5.14) 



with 



a = «f-z;f, + c = 1 + + (5.15) 



and that a general solution of (|5.14() is given by ^| eqs. (2.67)] 

Tj^^ = F(a,6,c;t) = (y -^^ '^u''~'={l ~ uY-'''\t - du. (5.16) 

Noting that in (|5?TT|l we have v[°^ + v^^^ = 0, it thus follows from ifTT^ . (|5?T^ . ifTT^ and (|5?T^ that 

f„^^ = det[,5J"+'=(t(t- l))^/2F(a,6,c;%fe=o,...„-i. (5.17) 
We know from ^17, proof of Prop. 6] that F satisfies the differential-difference relations 

t-^F{a,b,c;t) ^ -{c ~ b ~ l)F{a,b + I, c]t) - bF{a,b, c;t) (5.18) 
t{l-t)-^F{a,b,c;t) = {a - c)F{a - l,b, c;t) + {a - c + bt)F{a,b, c;t). (5.19) 

It follows that 

S{{t{t - l))''/^F{a, b, c; t)) = ^{t{t - Ijf^Fia, b,c;t) + {b + I - c)t'''\t - l)^'^+^F{a, 5 + 1, c; t) 
5({t{t - l)f/^F{a, b, c; t)) = (c - a - b/2){t{t - l)f/'^F{a, b, c; t) + (c - a){t{t - l)f'''^F{a - 1, b, c; t). 

^Froni these latter relations the working of the proof of |17[ Prop. 4] gives 

n-l 

fX^ = ]T(6+l-c),(c-a),t''"/2(<-l)''"/2+"(""i)/2det \F{a - 3,b + k,c;t)\ . (5.20) 

L J 7,/c— 0,...,n— 1 

J = l 

The method of the proof of 17, Prop. 5] allows this to be rewritten as a multiple integral, which when 
substituted in (|5.12l) and after substitution of the parameters according to (|5.15|) , H5.11|l shows 



'i-i pea r°° \ / r°° \ 



i=l l<j<k<n 

Replacing t by (is-l- l)/2, changing variables in the integrations uj i-^ {ivj + l)/2 and making use of H5.6|) 
we deduce that the proportionality constant in H5.10() (taken for convenience to be on the left hand side) 
is equal to 

JV-l 

(m n (1 + + ^),)*'^+'^^2-(''-^''-i)^7r"A%(,7,,y). (5.22) 
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Consequently, after recalling (|5.7|l . we conclude 



/ 1 vi( is + l 

V(1 + .52)A^M/2^JV { 2 



s = cot <j) / 2 , /.ti— >2/i 



7V-1 



( n (1 + * + + ^ + M),)^(^+^)^2-(2-^+i)^^^r^^[iV](e'^M, «^i,^i'2;r). (5.23) 



The use of this result lies with the fact that the r-function sequence {T]^'^(t; v)}Ar=o.i,... for general 
t and V = [yf'' , Uj'^'' , ^3°^ + N, v^f"^ ) satisfies, according to |17l proof of Prop. 2] , an equation of the form 

arm . 



= qN{l ~ qN)p% + 2vf\NPN ~ {vf~^ + I'f W + {^1^ + 4°^ + N){^T - vf^ + (5.24) 



and in addition recurrences determining {pn}n=i,2,... and {qn\n.=i,2.... are also known. Regarding the 
latter, set 



9n 



9n - 1' 



fn := g„(g„-l)p„ + (l + n-4°'-af )(g„-l)-4"^g„-(4")+n)*ii^^^ (5.25) 

where 

4") + =„f)-4°), af . (5.26) 

Then we have [13 Prop. 10] 

i (/„ + l + n- a^^'>)ifn + 1 + n - - 4°^) 



/n(/n + af) 

+ = -4°^ + ^^ + > , (5.27) 



which are a version of the discrete Painleve V equations |S0|. To use these recurrences to determine {pn}, 
{qn} given /o, go we first iterate (|5.27|) to determine {/„}, {(?«}. According to the first equation in (|5.25l) 
g„ can then be calculated in terms of gn, 

qn = (5.28) 

9n ~ 1 

Now that is known the second equation in H5.25|l allows p„ to be calculated in terms of /„, g„, 

Pn = -T^—T-, (fn + {I + n - 4"' - 4"')(1 - 9«) + 4"''?n + (4"' + uf-^^^^^) . (5.29) 
qn\qn ^} ^ qn 

To calculate {T^^[iV](e"^; /x, wi, i(j2; C)}jv=2,3.... the following recurrence scheme can therefore be given. 

Proposition 9. Let t^-'IN] T^"'^[-/V](e"^; /i, wi, W2; as specified by and letpM, qN denote the 

conjugate variables in the Hamiltonian i|5.ff)) to^/i parameters ^5.11]) . We have 

9n 



qn 



gn - 1 



^" 9n(gn - 1) 



- — rrf/n + (1 + ^^ + w + /i)(l - g„) - {w + ij)qn 

'L -L ) ^ 

g»(gr» - 1) 

qn - (1 - e*-^)- 
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where {fn}n=o,i,..., {5n}n=o,i,... flJ^e determined by the recurrences 

. (/„ + ! + n)(/„ + l + n + w + fi) 

9n+ign = e — — 

/«(/« - w - Aij) 

/„ + /„-! = + M + -— T - . _ (5-31) 

subject to the initial conditions 

50 = ^^, fo^{l + w + fi){qo ~l) + {w + fi)qo - {2w^ + 1) '^"^f " (5.32) 
go - 1 go - (1 - e"P) ^ 

Given t^^[0] = 1, and t^^[1] as t/ie element of the determinant in ^d.SJjl with j — k = Q, we have that 
{T^''[N]}pf=2,3,... can be computed in terms of {pn}n=i.2,... and {qN}N=i.2,... by the recurrence 

,t^'[N+1]t^^[N-1] 
{N + w + p)iN + w + f,) (^v-/[^])2 

= qN{l - qN)p% - 2{wi + fi)qNPN + {w + fi)pN + N{N + 2wi + 2fi). (5.34) 

Proof. The only remaining point to require explanation is the initial conditions (|5.32|) . H5.33|l . These 
come about because the PVI system admits the solution |EI Prop. 3] 

po = 0, t{t-l)j^\ogT];\t)^~ii{qo-t). (5.35) 

According to H5.23|l we require t = ^{is + 1), s = cot 0/2 and so < = 1/(1 — e*"^). Hence the second 
equation in H5.35(l . together with H5.23(l in the case iV = 1, gives (I5.33() . Also, setting n = in the second 
equation of H5.3()|l and equating the right hand side to zero gives the second initial condition in H5.32|l . 
The first initial condition in H5.32(l follows immediately from the first equation in l|5.30(l . □ 

^From the Okamoto theory 26 we know p — p{t; a), q — q(t; a) must satisfy a number of transfor- 
mation formulas with respect to t and a. Thus with 

:= (ao, ai, q;2, q;4, as) := (ao, 014, ^2, ^3, "i) := ("4, ai, "2, as, "o) 

one has 

p{l — t, a.^) — --p{t; a), q{l — t, a^) = 1 — q{t; a.) 

^(7; a^) = -a2q{t\ a) - q^{t] a.)p{t; a), g(y; a^) = / 
\t / t q\t\a.) 

a^) =-{t- l)p{t; a), g(^; = (5.36) 

Setting t = 1/(1 — e*"^) and inverting these formulas we could, if required, write down a variant of 
Proposition |51 in each of these cases which implicitly involves the variables 1 — t, l/t and t/{t — 1) 
respectively, but at the expense of permuting the a's. Consider in particular the first mapping involving 
t 1-^ 1 — t, a 1-^ a^. With t = 1/(1 — e**^) and as = —w — 11, = —w — fi this corresponds to simply 
taking the comlex conjugate. Indeed making the replacements p 1— *■ — P, g i— > 1 — Q we see that the right 
hand side of H5.34|l formally becomes equal to its complex conjugate, provided we identify P, Q with p, q 
respectively. 
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Let us now turn our attention to the special case 2/i e Z>o and ^ = 0. Then according to (|5.3|l . with 
t given by (|5.1|l i^^T^^[7V] is a polynomial in t, 



t'^^'T^'lN] = det 



r(2wi + 1) 

r{w - fi + j - k+ l)T{w + fi- j + k+l) 



X2Fi{-2fi,-w- fi + j -k;w- fi + j -k+l;t) . (5.37) 

J j,k=l,...,N 

Furthermore in this situation H5.33|l reduces to 

90 = 1 - — -77log2Fi(-2/i,-u;- /i + l;t) = — — (5-38) 

Zfj, at 2^i(— — w — /x; w — + 1; rj 

Also for 2/i <E Z>o we have the duality type relation between averages |17l eq. (3.42)] 

/Y[^l^^-n.)/2^^^^^\m+m(^l + tzifA oc ( J] ^'^^''^^^'ll + (1 + (1 - • (5-39) 

Recalling 1)1. 9|l we thus have 

T^^[iV](t; /i, wi, ^2; 0) cx T^^[2/*](l - t; N/2, 1 (u) + 2zi; + /i + TV); 0), (5.40) 

which is the PVI analogue of (|4.32|) and (|2.1t)|) . Setting t — Q shows the proportionality constant to be 
equal to 

(2/z)! Mjv(-/i + w,/J + w) 
TV! M2^i2wi+ N,-fi-w)' 

Each of the quantities qn,Pn, fn,gn in Proposition El can be written as a ratio of r- functions. Intro- 
ducing for convenience 

T'''[n]{t;fi,w,w;0 := t''>'T'''[N]{t; f,,wi,W2;0, (5-41) 

on the basis of exact tabulations with initial condition 1)5. 38|) we are led to the formulas 

f^'[n + l]{t; fi - 1/2, w + 1/2, w - 1/2; Q W(^; ^i,w,iu + 1; Q 

~ T^^[n + l]{t;n,w,w;C) f^^[n](i; /t - 1/2, u- + 1/2, tD + 1/2; ^ 

f^^[n+ l](t;/i,-u;,w;C) 



(5.42) 



Pn = 2/i(t - 1) 



/„ = -(n + l)i 



'r^^[n](i;/i,u;,w;^) 

f^^[n](t; /^ - 1/2, w + 1/2, w + 1/2; Q f^^[n - l](t; /x + 1/2, w + 1/2, w + 1/2; Q 
f^^ [n] (t ; /i, w, + 1 ; f^^ [n] (t ; /i, w + 1 , iD; e) 

•^^^[n+ l](i;/i- l/2,w- l/2,u) + l/2;0 T^^[n]{t; fi,w + l,w;£.) 



[n + 1] (t; /(, w; ^ [n] (t; /( - 1/2, w + 1/2, w + 1/2; ^ 

(5.44) 

1 f^^[n + 1] [t- /i - 1/2, w + 1/2, 15 - 1/2; C) T^^[n] (t; /i, w, w + 1; ,e) 



t fVi[n + l](i; /i - 1/2, w - 1/2, w + 1/2; ^ f ^^[n](t; /i, w + 1, iD; ^ ' 



(5.45) 



We have not completed a proof of these relations. However, as with the formulas (|4.33(l . an outline of 
how one goes about proving (|5.42|l - (|5.45|l is provided by the proof of Proposition 01 Here the matter is 
complicated by there being four fundamental Schlesinger operators instead of the two in PHI' theory, and 
the fact that p„, qn are functions of 1/(1 — t) rather than t. Some details of dealing with the first of these 
complications is given in |23| . while use of the transformation identities H5.36|l is required to deal with 
the second. Such arguing can only be used to establish 1)5. 42(1 - 1)5. 45|l up to proportionality constants. To 
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determine the latter, we proceed on the assumption that the proportionahty is independent of ^, allowing 
us to set i — 0- Then according to H5.37|l and the formulas of Proposition |51wc must have 

1 ^ ^ , 1^ 

Qn ~ 1, ~ -Qn, In ^ -(n+l), p„ ' 



' t-*0 t— ►oo t-^oo t^O W — II + n 

On the other hand it follows from (|5.41|) and the definition of Mi^{a, h) in H4.31|l that with ^ = 

Now using the evaluation formula in (|4.31|) we deduce the proportionalities in H5.42|I " H5.45|) . 

In relation to the PVI r-function f^^[rt] as specified by (|1.14l) we must first recall some theory from 
jl7| . Thus we know that the working leading to H5.21|l can be carried through with 

oo Jt ' ^Jo Jt 

and hence the PVI system admits a r-function sequence 

/ -i\ri(n-l)/2 



Comparison with (|1.14|l shows 

(^(« + b), \{b- a), -i(a + 6), -i(a + b)- /i)) 
^ ^_^^M(N-,y Ma^N^b-N) _ ^^^^^ _ , _ ^ _ ^. ^) (5.46) 



where 



JN{a,b) 



•^0 "^O l<j<k<N 

= l\ Tia + b+l + N + j) ^'-''^ 

It follows from this and (|5.24|) that we can compute {f^-' [N]{t; fi,a — N,b — N; (,)}n=q,i,... by a recurrence 
analogous to that in Proposition^ although we do not pursue the details. 



6 Applications 

In this section we will present results from the numerical evaluation of examples of the r-functions p.7|) 
and p.9|l based on the recurrences of Propositions |21 and IHl respectively. Consider first the r-function 
p.7|l . The cases = and jj, — 2 have particular significance. Thus let i?Ar_2(0, (0, s);a:"e~^) denote the 
probability that there are no eigenvalues in the interval (0, s) of the LUEjv as specified by the eigenvalue 
probability density function (|1.15|l . and let pN,2{0, {0, s); x°-e~^) denote the probability density of the 
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smallest eigenvalue in the same ensemble. These two quantities are inter-related by a single differentiation, 

PnAO, (0, s); a:'^e-") = --^Sjv,2(0, (0, s); ai^e""). (6.1) 

as 

To make contact with (II. 7|) consider the scaled limit of these quantities, 
^2^"'-^(0,(0,t)) := lim i;jv,2(0,(0,</4iV);x'^e-") 



N- 



(the reason for the superscripts "hard" is that the neighbourhood of the origin in the Laguerre ensemble 
is referred to as the hard edge; see e.g. QH). Now we know from J5] that 

E'r''iO,iO,t)) = e-*/4det[/,_,(\/i)],,fe^i,...,, = e^^/V"' [a](t; 0) (6.2) 

P'r'(0,t) = ie-*/Met[/2+,_.(V^)],,fc=i,...,a = ie-*/V"^'H(i;2) (6.3) 

where in both cases the second equality follows from (|3.2I) . Analogous to (|6.1|) we have 

p'r'*(o,i) = ~^i?r^(o,(o,t)). 

The large a limit of (|6.2() . (|6.3() is particularly interesting. Thus according to the Baik-Deift- Johansson 
theorem ^ (see [H] for a recent simplified proof) 

lim E^^"^{0, (0,0^ - 2a{a/2y/h)) = £;rf*(0, (s,oo)) (6.4) 

a — >oo 

where £'2°^'(0, (s,oo)) denotes the scaled probability of no eigenvalues in the neighbourhood of infinity, 
and similarly 

lim (2a(a/2)i/3s)p2(0, - 2a{a/2f/^s) = p^/*(0, s) 

a — *oo 

where p|°^*(0, s) denotes the scaled distribution of the largest eigenvalue. Let us then address the task of 
computing 

E'r\0, (0, - 2a{a/2f'^s)) =: g'^^^\a; s) (6.5) 

using H6.2|l and the recurrence scheme of Proposition |21 First it is clear that for large a and s of order 
unity the sequence 

{e-*/V"'[n](t;0) ^ ^^^Jn=o.u....a 

consists initially of numbers very small in magnitude. Hence it is necessary to work with high precision 
arithmetic throughout the calculation to ensure an accurate final result for the final member, which is 

equal to g^'^'^'^ {a; s) . This sequence in turn is calculated in terms of the sequence {Pn}n=o,i a-i a-s 

specified by the recurrence in Proposition 13 with t replaced by t/4. For the specific value s = 0.5 the 
results of Table ^a-re thereby obtained. 
The data fits well the extrapolation 

g''-'^(a;0.5)=go + 4k + - 
giving go = 0.990543 and thus from (|6.4|l predicting 

E^°^\0, (0.5, oo)) = 0.990543. (6.6) 

In fact i?|°^*(0, (s, oo)) is known in terms of a particular Painleve II transcendent q{s) High precision 
data by way of the values of £^|°*'(0, (0, oo)), q{0),q'{0) to 50 decimals have recently been given |29j . 
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a 


5r'''''''^(a;0.5) 


60 


0.991338737 


80 


0.991201326 


100 


0.991111203 


120 


0.991046762 


140 


0.990997995 


160 


0.990959574 



Table 1: Tabulation of g^'''''^ (a; 0.5) as specified by (jHSJ in the case s = 0.5. 

allowing for accurate determination of £"1°^* for gener al s. One finds £;|°f*(0, (0.5, oo)) = 0.990544..., 
showing us that Ht).6|) is accurate to 1 part in 10^. 

We now turn our attention to a particular example of the t- function H1.9I) . Let p'^J^ q{9) denote the 
probability density function for the spacing between consecutive eigenvalues in the CUE at or equivalently 
U{N). Then as noted in it follows from the definitions that 

.. . 2^X r^^[jV-2](e^"^/^;l,l,0;l) 
-[N -l)sm — ^v,[^_2](l;l,l,0;l) ' ^^'^^ 













^ N J 



Use of (|4.31|l shows 



.->](l;l,l,0;l)=(^^)^(^ + ^)(^ + ^) 



12 

and thus 

{^y--%^) - 4^ ^ -"[^ - 2Ke^-"/"; 1, M; 1). (6.8) 

According to Proposition |51 the key quantity in computing {T^^[ri]} by recurrence is t^^[1]. Now, from 
(ESI) and (ESJ) 

T^^(e^^l,l,0;l) 

= e-^^i^^i(-2, -2; 1; e^'^) + j^e-'^il - e'*)^Fi(3, 3; 6; 1 - e'*) (6.9) 

(identities for the 2F1 function can be used to check that this quantity is real. Using this in Proposition 
El and again using high precision computing, for the specific value X = 1/10 and a sequence of N values, 
we evaluated (16.8(1 . obtaining the data listed in Table 12 
The limiting distribution 

p^HX) = hm Ip^^-'I(^) (6.10) 

can itself be expressed in terms of a Painleve transcendent [201 ^] • Moreover its power series about 
X = is known to high accuracy ||19| , and from this we can compute 

pbuik/j^) = 0.032468767196387... (6.11) 

Extrapolating the data of Table [5] using the ansatz 

27r Q-^^ 2nX si S2 

gives So = 0.032468767193... which agrees with l|HTT|l to 3 parts in 10^^ 
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AT 

iv 


27r„CUE^27rXN, 
nPN-2\ N ) 


X=l/10 


10 


0.03215040321 


30 


0.03243339939 


50 


0.03245603495 


70 


0.03246227118 


90 


0.03246483751 



Table 2: Tabulation of the scaled probability density at X = 0.1 for the spacing between 
consecutive eigenvalues in the CUEiv- 
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